Heat flow in the postquasistatic approximation 



O 

< 
in 



o 

cr 
I 

$-1 



> 

in 

o 
o 



X 



B. Rodriguez-Mueller,^ C. Peralta,^''^ W. Barreto,'* and L. Resales^ 

^ Computational Science Research Center, College of Sciences, 

San Diego State University, San Diego, California, USA 

^ Deutscher Wetterdienst, Frankfurter Str. 135, 63067 Offenbach, Germany 

^School of Physics, University of Melbourne, Parkvtlle, VIC 3010, Australia 

''Centro de Fisica Fundamental, Facultad de Cienctas, Universidad de Los Andes, Merida, Venezuela* 

^ Laboratono de Fisica Computacional, Universidad Experimental 

Politecnica "Antonio Jose de Sucre", Puerto Ordaz, Venezuela 

(Dated: August 6, 2010) 

We apply the postquasistatic approximation to study the evolution of spherically symmetric fluid 
distributions undergoing dissipation in the form of radial heat flow. For a model which corresponds 
to an incompressible fluid departing from the static equilibrium, it is not possible to go far from 
the initial state after the emission of a small amount of energy. Initially collapsing distributions 
of matter are not permitted. Emission of energy can be considered as a mechanism to avoid the 
collapse. If the distribution collapses initially and emits one hundredth of the initial mass only the 
outermost layers evolve. For a model which corresponds to a highly compressed Fermi gas, only the 
outermost shell can evolve with a shorter hydrodynamic time scale. 

PACS numbers: 04.25.-g,04.25.D-,0.40.-b 



I. INTRODUCTION 

Dissipation as an emission process [T] is crucial for the 
outcome of gravitational collapse. Thermal conduction is 
usually considered proportional to the gradient of tem- 
perature. This is a sensible choice, since the mean free 
path of particles responsible for the propagation of en- 
ergy in stellar interiors is very small as compared with 
the typical length of the object [2|. Observations from 
supernova 1987A indicate that the regime of radiation 
transport prevailing during the emission process is closer 
to the diffusion approximation than to the free stream- 
ing limit [3 . The addition of a test bed for studying 
dissipation mechanisms and other transport processes in 
order to later incorporate them into a more sophisticated 
numerical framework (Arnowitt-Deser-Misner [ADM] or 
characteristic) is a necessity. 

In this work we study a self-gravitating spherical dis- 
tribution of matter containing a dissipative fluid in the 
diffusion limit. We found behaviors similar to those re- 
ported with a different mechanism by [Ij , and report the 
zeroth order results for dissipation. We use noncomov- 
ing coordinates and follow the method reported in |S] , [5] 
named the postquasistatic approximation (PQSA), which 
has been proposed as a test bed in numerical relativity 
[7]. For recent advances and applications see and [9]. 
For origin, reviews and details of the PQSA see [10] , [TT]- 
[T5] and |Ilj-[in]- We do not consider here temperature 
profiles to determine which processes can take place dur- 
ing the collapse. For that purpose, transport equations 
in the relaxation time approximation have been proposed 
to avoid pathological behaviors (see for instance [5D] and 



'On sabbatical leave. 



references therein). These issues will be considered in 
a future investigation. It is worth mentioning here that 
in order to get a higher order approximation we have 
to know the zero order approximation in the relaxation 
time, as in the present study. To the best of our knowl- 
edge, no author has undertaken in practice the dissipative 
matter problem in numerical relativity. Our purpose here 
is to show how heat flow processes can be considered in 
the context of the PQSA. The results indicate that an 
observer using radiation coordinates does not "see" some 
details when heat flow is considered. The final goal is 
to eventually study the same problem using the Miiller- 
Israel-Stewart theory for the dissipative system, which is 
highly nontrivial in spherical symmetry. 

In Sec. II, we present the field equations, the match- 
ing conditions and the set of surface equations. For ad- 
ditional details concerning the PQSA method see [^ and 
[3]. Three models are presented in Sec. Ill and some 
remarks are discussed in Sec. IV. 



II. MAIN EQUATIONS 

To write the Einstein field equations we use the line 
element in Schwarzschild-like coordinates 

ds^ = e^dt^ - e^dr^ - r^ {dd^ + sin ^ed(f>^) , (1) 

where v — t^(t, r) and A = A(i,r), with {t,r,9,<j)) = 
(0,1,2,3). 

In order to get physical input we introduce the 
Minkowski coordinates (r, x, y, z) by [H] 

dr = e^/^dt, dx = e^/'^dr, dy = rd0, dz = r sin 6*^0, (2) 

In these expressions i^ and A are constants, because they 
have only local values. 



Following the Bondian point of view as in [5T] , [B] , [7] 
and [4] we assume that, for an observer moving relative 
to the local Minkowskian coordinates with velocity w in 
the radial direction, the space contains an isotropic fluid 
of energy density p, radial pressure p, and radial heat 
flux q. For this comoving observer, the covariant energy 
tensor in Minkowski coordinates is thus 
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Making a Lorentz boost we write the field equations in 
relativistic units (G = c = f ) as follows [4]: 



f 

87rr 



87rr 



f 



f 



A. 



e " \ - + V r 
r 



(4) 



(5) 



1 2 

62tt ' r 



e~''[2X,tt + X,t{\t~'^,t)]}^ 



S = 



87rr 



(6) 



(7) 



where the comma (,) represents partial differentiation 
with respect to the indicated coordinate and the con- 
servative variables are 
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and the flux variable 
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as in the standard ADM 3+1 formulation. Within the 
PQSA p and p are referred as to effective density and 
effective pressure, respectively. Note that from (pi) the 
velocity of matter in Schwarzschild coordinates is 



dr 
di 



= c.e('^-^)/^ 



It is easy to check that |2] 

Pa = Qa, 



(11) 



(12) 



which expresses the continuity of the radial pressure 
across the boundary of the distribution r = a{t). Equiv- 
alently, in terms of the effective variables 



PaW^ 



Pa = Pa^a+lai^ 

Defining the mass function as 

e"^ = 1 - 2m/r, 
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(13) 



(14) 



and substituting (14 1 into (|4| and ^ we obtain, after 
some rearrangements. 
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This equation is the momentum constraint in the ADM 
3+1 formulation, it expresses the power across any mov- 
ing spherical shell. 

Equation (|6| can be written as T^ = or equivalcntly, 
after a lenghiy calculation 
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This last equation is the generalization of the Tolnian- 
Oppenheimer-Volkov for nonstatic radiative situations. 
In can be shown that Eq. ( 16 ) is equivalent to the equa- 



tion of motion for the fluid in conservative form in the 
standard ADM 3+1 formulation 17 . 



At the surface, Eqs. ( 11 ), ( 15 ) and ( 16 1 lead us to a set 



of differential equations for a, ma and uja if we prescribe 
in some way the metric functions (m, and v). 

The other two field equations (|4| and ^ can be inte- 
grated to obtain 



TO = / 47rr p dr 
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which is the Hamiltonian constraint in the ADM 3+1 
formulation and 
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the polar slicing condition, from where it is obvious that 
for a given radial dependence of the effective variables, 
the radial dependence of the metric functions becomes 
completely determined. 



III. MODELING 

We consider here a seed model inspired by the well- 
known Schwarzschild interior solution. This model cor- 
responds to an incompressible fluid departing from the 
static equilibrium. Following the PQSA we take 
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FIG. 1: Evolution of the radial local velocity uj (multiplied by 
10'^) for the Schwarzschild-like model. The initial conditions 
are a(0) = 5.0, m(0) = 1.0, uJa{0) = 0.0. The total radiated 
mass is Mr — 10^*ma(0). 



FIG. 3: Evolution of the radius a for the Schwarzschild- 
like model. The initial conditions are a(0) — 5.0, m(0) = 
1.0, cja(O) = -0.001. The total radiated mass is Mr = 
10-^771,(0). 
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FIG. 2: Evolution of the heat flow q (multiplied by 10^) for the 
Schwarzschild-like model. The initial conditions are a(0) = 
5.0, 771(0) = 1.0, uja{0) ~ 0.0. The total radiated mass is 

Mr = lO""* 771,(0). 



1.917 



1. 915 



^ 1.913 



1. 911 



1.909 




0.0 2.5 5.0 7.5 10.0 

Time 

FIG. 4: Evolution of the energy density p (multiplied by 10^) 
for the Schwarzschild-like model. The initial conditions are 
a(0) = 5.0, m{0) = 1.0, cJa(O) = -0.001. The total radiated 
mass is Mr = 10"^ 777,(0). 



where / is an arbitrary function of t. The expression for 
p is 
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where fc is a function of t to be defined from the boundary 



condition ([12|) or (jlj). Thus, (|20|) and (|13|) give 
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with 
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FIG. 5: Evolution of the radius a for the Tolman Vl-hke 
model. The initial conditions are a(0) — 8.0, m(0) = 1.0, 
c»Ja(0) = -0.02. The total radiated mass is Mr = 10"" ma{0), 
with a narrow Gaussian given by E = 0.01 with maximum at 
to = 1.0. 
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Using (17) and (18) it is easy to obtain expressions for 
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Thus, the system of equations at the surface can be in- 
tegrated, but it is necessary to specify one function of t 
and the initial data. We choose 



L = 4na' 
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(27) 



We explore a complete range of initial conditions and 
parameters of integration to get physically acceptable re- 
sults. A radiated mass bigger than 10~^ and an initially 
contracting velocity are not permitted. The reason is 
a complex root calculating the local radial velocity and 
heat flow for some regions of spacetime. Representa- 
tive and acceptable results are shown in Figs. 1-2. For 
this model, the energy density and the radius of the dis- 
tribution remain almost constant (within six significant 
figures). These features were not reported in the past 
using radiation coordinates and lead us to the follow- 
ing model. Our results clearly show that the heat flow 
keeps the evolution near quasistaticity (slow evolution). 
Under the same initial compactness used above, that is, 
a(0) = 5, we found a possible initial local radial velocity 
of Wa(0) = -10"'^ and a radiated mass of Mr = 10"^, 
producing now an appreciable change in the energy den- 
sity and the radius of the distribution. These results 
are shown in Figs. 3-4. For these conditions only the 
evolution of a bubble is possible [r/a ~ 0.99 — > 1.00). 
We do not observe any evidence of thermal peeling [2], 
that is, positive velocities (expansion) of outer shells and 
negative velocities (contraction) of the inner shells. The 
development of thermal peeling leads to complex roots 
for the radial velocity. 

We consider now other interior seed model based on the 
Tolman VI interior solution _22] . This model corresponds 
to a highly compressed Fermi gas. Let us take 
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where g and a are functions of i, which can be determined 
using (13). Thus 
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with Lo = M,./VS7r, to = 5.0 and E = 0.25, which corre- 
sponds to a pulse radiating away a fraction of the initial 
mass Mr- Therefore, the system can be numerically in- 
tegrated for the following typical initial conditions: 

a(0) = 5.0, TOa(O) = 1.0, UJa{Q) = 0-0. 

The integration was done up to some t guaranteeing well 
behavior of the physical variables, that is, p > 0; p > p; 
1^1 < 1; w, g G 5ft. Feeding back the numerical values 
of a, ma and w^ (and their derivatives) in (17) and (18) 
we obtain m and v (and their partial derivatives) for any 
value of r. Thus, variables p, p, w and q can be monitored 
for any piece of the material, via field equations. We cal- 
culated them for the values r/a = 0.0, 0.2, 0.4, 0.6, 0.8 
and 1.0. 
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Once the metric functions are obtained from 
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17) and 



( 18 ), the system of equations at the surface can be agam 
numerically integrated for the following initial conditions: 

a(0) = 8.0, TOa(O) = 1.0, uja{0) = -0.02. 

As before, a radiated mass bigger than 10^'* is not per- 
mitted. But even more, now it is not possible to go inside 
the distribution without violating real values assumption 
from the beginning. For that reason the Gaussian has 
been set to E = 0.01 and 1^ = 1. At the surface, see Fig. 
5, the results are as expected. 



IV. CONCLUSIONS 

In this paper we considered heat flow as a transport 
mechanism in the PQSA. Heat flow produces a stable 
configuration, which is the opposite effect of viscosity 
[3] . This result indicates that a combination of viscosity 
(anisotropy) with heat flow may be crucial for gravita- 
tional collapse or at least just out of equilibrium, where 
we expect the PQSA is a good approach. We did addi- 
tional tests including anisotropy but its effect is marginal. 
For distributions far from equilibrium we find that heat 
flow is a very restrictive transport mechanism. 

These results are apparently different for the same 
configurations in radiation coordinates [in]-[Il]. If the 
initial distribution is in equilibrium the transition from 
static to postquasistatic, in radiation coordinates, al- 



lows the sphere to "instantaneously" bypass diffusion 
stressing. But in Schwarzschild coordinates we can fol- 
low the transition from the static configuration to the 
postquasistatic with more resolution. Diffunding radi- 
ation strongly interacts with matter. As a result, the 
interior of the distribution is not permitted to go far 
from equilibrium. When including heat fiow, there is 
not PQSA except very close to or at the surface. 
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